Numerical solutions of the three-dimensional magnetohydrodynamic alpha-model 
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We present direct numerical simulations and a-model simulations of four familiar three-dimensional magne- 
tohydrodynamic (MHD) turbulence effects: selective decay, dynamic alignment, inverse cascade of magnetic 
helicity, and the helical dynamo effect. The MHD Q-model is shown to capture the long-wavelength spectra in 
all these problems, allowing for a significant reduction of computer time and memory at the same kinetic and 
magnetic Reynolds numbers. In the helical dynamo, not only does the a-model correctly reproduce the growth 
rate of magnetic energy during the kinematic regime, but it also captures the nonlinear saturation level and the 
late generation of a large scale magnetic field by the helical turbulence. 

PACS numbers: 47.27.Eq; 47.27.Gs; 47.11.+j 



I. INTRODUCTION 

The "alpha model," as it has come to be called in fluid me- 
chanics, is a procedure whereby, by suppressing small spa- 
tial scales in a computation in a way that intends to do mini- 
mum damage to the accuracy with which the long wavelength 
spectral components are calculated, one can realize substantial 
savings in computing tiniejni|l|l|i|l|S|l|i|i[i3[iil[i3 
fl3ll . In a previous paper Il4ll. we gave a simple way to extend 
the alpha model to magnetohydrodynamics (see also 0, Holl 
for extensions in the non-dissipative case), we specialized it to 
two dimensions, and numerically tested its predictions in a se- 
ries of computations. These were chosen as situations where 
direct numerical simulations (DNS) that started from identi- 
cal initial conditions were feasible. The intent of this present 
paper is to present comparisons of the same kind for three- 
dimensional (3D) magnetohydrodynamics (MHD). This is a 
straightforward program in the light of the two-dimensional 
(2D) investigations and we will draw heavily on the ma- 
terial in that paper to avoid repetition, but in 3D, new phenom- 
ena arise, such as the generation of magnetic fields through 
stretching by velocity gradients, and furthermore 3D is com- 
putationally more demanding than 2D. 

In Sectionini we take the 3D alpha model MHD equations 
fl4ll and describe briefly four problems upon which they will 
be tested against DNS treatments of the same problems. They 
are selective decay, dynamic alignment, the inverse cascade of 
magnetic helicity, and the mechanically driven turbulent dy- 
namo. The first two have already been tested in 2D fTl l and 
the third has a 2D analogue in the inverse cascade of magnetic 
vector potential kMI . The fourth is also not an unfamiliar ef- 
fect, and we have recently been involved in addressing it for 
the special case of low magnetic Prandtl number ITsIl and for 
non-helical flows. 

Our conclusions reached in Sees. IllllVII are consistent for 
the most part with those reached for 2D MHD: the alpha 
model does a good job of reproducing the spectral behav- 
ior of the long-wavelength Fourier amplitudes (wavenumber 
k < a^^, where a is the spatial scale over which the velocity 
field and magnetic field are filtered). Because of the relative 
lack of surprises in the selective decay and dynamic alignment 
Sections, we rely on relatively brief presentations, to then fo- 



cus in the study of the inverse cascade of magnetic helicity 
and the dynamo effect. Finally, we summarize the results in 
Sec.rvlll 



II. RELEVANT EQUATIONS; PROBLEMS CONSIDERED 

In familiar "Alfvenic" dimensionless units, the original 
MHD equations are 



dv 

at 

OB 

'dt ' 



-v-Vv = -VP +j X B - j/V X (1) 
V • VB B ■ Vv - r;V X j, (2) 



together with V • v = = V • B. 

The velocity field is v, the magnetic field is B = V x A, 
where A is the vector potential. The electric current density is 
j = V x B and the vorticity is a; = V x v. The dimensionless 
pressure, normalized to the (uniform) mass density is V, and 
is obtained by taking the divergence of Eq. O, using the in- 
compressibility condition V ■ v = 0, and solving the resulting 
Poisson equation. Removing a curl from Eq. (|3 gives 



dA 

— - = V x B - ?;j - V$, 
ot 



(3) 



where <I> is the scalar potential, obtainable also from a Pois- 
son equation by imposing the Coulomb gauge V • A = 0. 
The kinematic viscosity is and the magnetic diffusivity is 
f]. In these dimensionless units, can be interpreted as 
a Reynolds number = UL/v where in laboratory (c.g.s) 
units, [/ is a mean flow speed and L is a length character- 
istic of it. Similarly, 77^^ can be interpreted as a magnetic 
Reynolds number i?„i = UL/rj. The magnetic Prandtl num- 
ber is P,n =v/l] = Rm/ Re- 
in the alpha model the fields v and B are smoothed but a> 
and j are not iHIT^l . The prescription is 



B, 



^3^, exp[-|x-x-|/a] ^ (4) 

iTTa^lx — x'l 
^3^, exp[-|x-x-|/a] ^ (5) 

47ra^ X — x' 
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Here a is an arbitrary filtering length, generally chosen 
smaller than the length scales one wishes to resolve. If v and 
B are Fourier-decomposed 



(6) 
(7) 



the connection between the Fourier transforms of the 

smoothed fields Ug and and Vk(t), Bk(t) are 

u,(k,i) = Vk(t)/(l + fc2a2) 

B«(k,0 = Bk(t)/(l + fc'a'), (9) 



or in configuration space 

B = (l-a^V^)B,. 



(10) 
(11) 



Note that we choose to smooth both the velocity and the mag- 
netic field at the same length-scale, a choice appropriate for 
the unit magnetic Prandtl number (j/ = 77) cases treated in this 
paper (for a different choice, see llSll '). 

The dynamics of the alpha model fl4ll amount to solving 
the pair, 

^ + u, • Vv = -vjVui -VP+jxB, 

-vV X (12) 

9Bc 

— ^ + U3-VB,, = B, -Vu, -TyVxj, (13) 
at 



where it is to be emphasized that in Eqs. (I12I13> . v, j, and uj 
are not smoothed. V is to be determined, as before, from the 
relevant Poisson equation. 

In rectangular periodic boundary conditions (which we em- 
ploy throughout), the ideal (v ~ Q ~ if) invariants that have 
been identified for Eqs. (I12I13> are the energy E 



v + B -B^d^x, 



the cross helicity He, 

He 



- I v Bsd-'x, 



and the magnetic helicity Hm 
1 

Hm = - 



A, • B, d^ 



(14) 



(15) 



(16) 



In the presence of non-zero rj and i/, the decay rates for E, 
He, and Hm can readily been shown to be 



dE 

'dt 
dHc 

dt 
dHM 

dt 



= —V 
1 



ujs-i^d^x-'q I fd^x (17) 



-u / uj-isd'x~ -ij I u jd^x (18) 



When we write or A,,, we mean that the same smooth- 
ing recipe as in Eqs. j4l5> has been applied to the unsmoothed 
fields u} or A. It is possible, and sometimes desirable, to use 
different smoothing lengths a^, as for the mechanical and 
magnetic quantities 1 14. isil . 

E, He and Hm as defined here are the ideal invariants of 
Eqs. ( I12I13> . and reduce, as a ^ 0, to the usual ideal 3D 
MHD invariants. Sometimes, to make the global quantities 
agree at i = for initial-value runs, we may initially normal- 
ize the Fourier coefficients of the fields, by multiplication by 
a common factor, to bring E and He into exact initial agree- 
ment with the corresponding numbers for the ideal 3D MHD 
invariants (note that Hm involves two smoothed fields, and 
therefore can not be matched to the DNS initial conditions at 
the same time). Hereafter, such global quantities as E, He, 
etc., will be referred to unit volume. 

It is well known that for decaying turbulent situations, the 
presence of enough initial Hm or He can lead to a late- 
time state in which the ratios \Hm/E\ or \He/E\ can be 
close to maximal. The first situation, called "selective decay," 
leads to a late-time quasi-steady state in which 
the remaining energy is nearly all magnetic and is nearly all 
condensed into the longest wavelength modes allowed by the 
boundary conditions. The second situation, called "dynamic 
alignment leads to a late-time quasi-steady state 

in which v and B are nearly parallel or anti-parallel. In both 
cases, the states can be very long-lived because the nonlinear 
transfer to small scales has essentially been shut down ("sup- 
pression of nonlinearity"). We illustrate these two situations 
in Sees. imiandHVl 



Inverse cascade processes l23il24l25Ll2ffl are those wherein 
excitations externally injected at the small scales are preferen- 
tially transferred to the larger scales and pile up there, creating 
coherent macroscopic structures at large scales where none 
were present initially. A quantity which can be inversely cas- 
caded in 3D MHD is Hm fZ^U^- We illustrate this with an 
externally-driven run in Sec. |V] 

Dynamo processes (see Ref. IztIi for a review) are those 
whereby mechanical injection of excitations transfer energy 
to magnetic fields, causing them to amplify. A novel example 
of helical dynamo action using the alpha-model is treated in 
Sec.|VT| 

In all four cases, well-resolved DNS solutions are regarded 
as baseline truths against which alpha-model computations are 
to be tested. 



= / j-Bsd^x 



(19) 



III. SELECTIVE DECAY 



In selective decay s, energy de cays rapidly relative to mag- 
netic helicity, if any I17l[l§.[l^ . In order to display the pro- 
cess most clearly, it helps to start an initial-value decay run 
with a significant amount of magnetic helicity. One way to 
accomplish this is to make the initial values of v and B out of 
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what are called "ABC" flows. We define 



'Vabc = [B cos(fcy) + C sin(A:2;)] x + 
+ [A sin(fca:;) + C cos{kz)] y - 
+ [Acos{kx) + Bsm{ky)] z 



(20) 



for arbitrary real numbers A, B, C, and k. vabc is an eigen- 
function of the curl. The specific initial conditions chosen are 

v(i = 0)= vo[^ABc{k,(l)k)+Hk)e'^-''] (21) 

k—kbot 

v^Bc(fc,'/'fc) + B(k)e^H . (22) 



B(t = 0)- ^ bo 



The notation "VABcik, <pk) means that for each k in the sum- 
mation, a random phase <f>k is added to the arguments of the 
sines and cosines for that k. The summations are over all the 
k values (which lie on a lattice in k space defined by the pe- 
riodic boundary conditions) between radii kbot and ktop- The 
v(k) and B(k) represent added random perturbations. 

Here, we have chosen A = B = C = 1, fcf,ot = 6, ktop ~ 
10, and vq, bo are chosen to make the initial ^w^) = (-8^) = 
1, where "(.)" means a spatial average over the basic box. It is 
also the case that initially, (v • B) = 0. Random modes v(k) 
and B(k) are added with an energetic level to initially give 
(A • B) = 0.5 (|A||B|). The dimensionless inverse Reynolds 
numbers are v ~ rj ^ 0.002. 

Three runs for a typical case are displayed. The first of these 
is a well-resolved DNS run at a resolution of 256^, with de- 
aliasing achieved by zeroing out all Fourier coefficients with 
k > 256/3, a method that will be used throughout (usually 
referred to as the "2/3 rule"). Then two a-model runs are 
performed with the same initial conditions, a 128^ run with 
a ~ 1/20 and a 64'^ run with a = 1/10. The same values 
of i/, 1] apply to all three runs. The caption of Fig. iden- 
tifies the decaying energies (kinetic energy Ek and magnetic 
energy Em) as functions of time. Fig. [^b shows the ratio 
(A ■ B) / (|A| |B|) as a function of time for the three runs; it 
has increased to above 0.999 by the final time. 

Fig. |5]shows the (unnormalized) energies and magnetic he- 
licities for the three runs. Note that by normalizing the DNS 
and a-model initial conditions to have equal energies, it has 
meant that the a-model magnetic helicities have necessarily 
started at lower initial values than those of the DNS. 

Figs. |3and|?]show the associated energy spectra plotted vs. 
wave number. Fig. |3]is at an early time {t ~ 10) and shows the 
total energy spectrum compensated by Kolmogorov's —5/3 
law. Fig. |4]shows kinetic (Ek) and magnetic (Ej^j) energy 
spectra at a very late time [t = 733). The two values of a^^ 
are shown as vertical lines. Below k ^ a^^, the DNS and 
a-model agree reasonably well. 

As follows from Figs. |2]and0] at late times the magnetic 
field is concentrated at large scales (k = 1) and has maximum 
relative helicity (note that E ^ Em ^ Hj^i after t ~ 200 
in both the DNS and alpha-model simulations). Fig. |5]shows 
surfaces of constant Hm ait — 800 in the 3D domain, for the 
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FIG. 1: (color online) (a) Magnetic energy (upper curves) and kinetic 
energy (lower blue curves) as a function of time until t = 5, and 
(b) relative magnetic helicity as a function of time until t = 600, 
for the selective decay runs. Solid lines correspond to DNS, dashed 
lines to 128'' a-model simulations, and dotted lines to 64'' a-model 
simulations. 



DNS and the 64^ alpha-model simulation. The alpha-model 
is able to reproduce the large scale structures observed in the 
DNS, and only slight differences can be observed. As will 
be shown in Section HVl this is not always the case when using 
periodic boundary conditions (similar results were obtained in 
2D MHD simulations fll v 



IV. DYNAMIC ALIGNMENT 



In this case, we load Fourier coefficients into the spherical 
shells with fc^ot = 6 < /c < ktop = 10 with equal ampli- 
tudes but enough correlation between the phases of v and B 
so that initially (v • B) ~ 0.3 (|v||B|); otherwise the phases 
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FIG. 2: (color online) Total energy (upper blue curves) and magnetic 
hielicity (lower curves) as a function of time. Labels are as in Fig.Q 
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FIG. 3: Total energy spectrum compensated by Kolmogorov's —5/3 
law, for the three dynamic alignment runs (labels are as in Fig. 0, 
att — 10. Vertical dotted and dashed lines indicate respectively the 
scales a^^ = 10 and 20. 



are random. We again do a 256'^ DNS run, an a-model run 
at 128'^ with a — 1/20, and another a-model run at 64'^ with 
a = 1/10. For all three runs, v = rj = 0.002. The same 
conventions are adopted for the graphics as in Sec. |lll] 

Figs. |6li,b show the decay of the kinetic and magnetic en- 
ergies (a), chosen initially to be equal; and (b) the degree 
of alignment, as measured by the mean cosine of the align- 
ment angle, (u,, ■ B) / (|us||B|) that develops as a function 
of time. Since much of the alignment is contributed by the 
small scales, the a-model underestimates the degree of align- 
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FIG. 4: (a) Kinetic and (b) magnetic energy spectra, for the three 
dynamic alignment runs (labels are as in Fig.0, att = 733. 

ment, and the disparity becomes greater as a^^ is decreased, 
though the accuracy remains within the 10 percent level. Fig. 
shows the decay of both E and He, with the more rapid 
decay of the former. There is, in this case, no preferential 
migration of any global quantity to long wavelengths. 

Figures|8]and|9]show the kinetic and magnetic energy spec- 
tra at an early time, t ~ 4.5, and at a late one, t = 156. The 
agreement of the a-model and DNS for k < a^^ is again seen 
to be excellent except for an unexplained over-estimate at the 
earlier time t — 4.5 for the kinetic energy spectrum. 

Fig. no| shows surfaces of constant He at t = 150 in the 
3D domain, for the DNS and the 64'^ alpha-model simulation. 
While there are marked similarities in the kinds of structures 
present in the DNS and in the alpha runs, there are no one- 
to-one correspondences as to specific features, either as to lo- 
cation or orientation. As in the 2D case fl4il we conclude 
that in this case the alpha-model does an excellent job repro- 
ducing the statistical properties of the large-scale spectra, but 
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FIG. 5: (color online) Surfaces of constant magnetic helicity density 
att = 800 at 90% of its maximum value, for the DNS (above), and 
the 64^ alpha-model simulation (below). 
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FIG. 6: (color online) (a) Magnetic energy (upper curves) and kinetic 
energy (lower blue curves) as a function of time until t = 5, and 
(b) relative cross helicity as a function of time until t = 160, for 
the dynamic alignment runs. Solid lines correspond to DNS, dashed 
lines to 128'' a-model simulations, and dotted lines to 64'' a-model 
simulations. 



V. INVERSE CASCADES 



small-scale detailed phase information (such as the location of 
structures) is lost. 

The reason for the striking agreement bet'ween the a-model 
and DNS exhibited in Fig. |5| as contrasted with the disagree- 
ment shown later in Fig. ^|is that in the case of selective 
decay, both computations have found the same final state: the 
isotropic, maximum-helicity, k — 1 state. This state is the 
"ABC flow" with A, B, and C all equal. 



Inverse cascades of magnetic helicity, driven mechanically 
at the small scales, have long been known to be an efficient 
dynamo mechanism for generating large-scale magnetic fields 
I2al . Here, we try a different approach: we drive the magnetic 
field directly at small scales with a random forcing function 
that is a superposition of "ABC" flows [see Eq. j20H between 
wavenumbers fc = 8 and 9, and with no driving for the ve- 
locity. The randomness is introduced by randomly changing 
the phases in the trigonometric arguments of each ABC com- 
ponent with a correlation time of At = 1.25 x 10~^ (in all 
the simulations we discuss in this section, the time step is 
2.5 X 10^"^). A tiny seed velocity field is amplified some- 
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FIG. 7: (color online) Total energy (upper blue curves) and cross 
helicity (lower curves) as a function of time. Labels are as in Fig.|6| 



what, but the kinetic energy always remains well below the 
level of the magnetic energy throughout. We again exhibit the 
results of a 256'^ DNS computation, and a-model computa- 
tions with a = 1/20 and 1/10, with resolutions of 128'^ and 
64'^ respectively. 

Figs. ^2 show the time histories of the energies (a) and 
magnetic helicities (b) for the three runs. The rather abrupt 
phase jumps in the ABC flows give the lines a jagged appear- 
ance and it is sometimes difficult to identify which of the three 
runs is which. Suffice it to say that the two a-model runs ex- 
hibit the same features as the DNS runs, but with a time lag 
that is greater for the larger a. This is visible more clearly 
in Fig. where the magnetic helicity spectra for the three 
simulations are plotted as functions of k. The curves are the 
helicity spectra as functions of time. The lower levels of ex- 
citation are associated with earlier times. The times exhibited 
range from t = 30 to t = 72.5. The peak, once established, 
moves to the left with nearly the same speed in each case. The 
suppression of small scales, where the unsmoothed ABC flow 
is also unstable, may be responsible for the time lag. This 
time-lagged behavior is reminiscent of what happened in two 
dimensions with the inverse cascade of mean square vector 
potential j3- However, note that in three dimensions once 
the inverse cascade has been established, the growth rate of 
magnetic helicity is well captured by the alpha-model (Fig. 
II lib), indicative of a more local cascade (in scale). The power 
laws present in the s pect ra of magnetic helicity, and kinetic 
and magnetic energy 12511 are also well captured by the alpha- 
model. 
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FIG. 8: (a) Kinetic and (b) magnetic energy spectra, for the three 
dynamic alignment runs (labels are as in Fig.|SJ, at t = 4.5. Vertical 
dotted and dashed lines indicate respectively the scales oT^ = 10 
and 20. 



VI. THE DYNAMO 

The mechanically-driven dynamo, in which injected me- 
chanical energy is converted to magnetic energ y a t large 
scales, has long been a recurrent problem in MHD 12711 . Here 
we are able to show that the alpha model yields the same re- 
sults within acceptable accuracy as those of a DNS of the same 
situation (see Ref. IT^ for another case of recent interest). 

We begin with a velocity field which is again forced ex- 
ternally with the "ABC" geometry of Eq. ( I20l i. We choose 
A = 0.9, B =. 1.0, and C = 1.1, /sq = 3, with ^ = 
V = 0.002. This choice is governed by the knowledge that 
the A = B — C flow gives the largest dynamo growth rate 
|2^ but it is hydrodynamically very stable Isoll : breaking that 
symmetry allows for turbulence to develop faster f3l|. 

The force is allowed to operate until a statistically-steady 



7 




10 
k 



FIG. 9: (a) Kinetic and (b) magnetic energy spectra, for the tliree 
dynamic alignment runs. 

turbulent Navier-Stokes flow prevails. Then a magnetic seed 
field is introduced at a very low level in the modes from k ~ 1 
to fc = 10. As in some other sections, we compare a DNS 
run at resolution 256'^ with two a-model runs, one with a = 
1/20 and 128^ resolution, the other with a = 1/10 and 64^ 
resolution. 

Before embarking on the MHD comparison between DNS 
and alpha-model results, it is instructive to compare the hydro- 
dynamic properties of the flow. When the small magnetic seed 
is introduced, the Lorentz force in the Navier-Stokes equation 
can be neglected. The induction equation is hnear in the mag- 
netic field, and as a result, the geometrical properties of the 
flow are responsible for the observed amplification. 

The flow generated by the external ABC force is helical. 
Previous studies of the alpha-model behavior in simulations 
of hydrodynamic flows were carried for non-helical flows 
0, flUl . As a result, here we will focus only on the charac- 
terization of the flow helicity. The amount of helicity in a 




FIG. 10: (color online) Surfaces of constant cross helicity density at 
t — 150 at 50% of its maximum value, for the DNS (above), and the 
64'' alpha-model simulation (below). 



flow (both for DNS and the alpha-model |7J|) is measured by 
the kinetic helicity 

HK = ^Jyru;d''x. (23) 

It is also useful to normalize this quantity introducing the rela- 
tive helicity 2Hk /{{\v\) (IwD). Fig.ll3lshows the probability 
distribution function (pdf) of relative kinetic helicity for the 
DNS and alpha-model simulations. A stronger positive tail 
can be identified in all cases, giving rise to a net positive ki- 
netic helicity in the flow. 

In 3D hydrodynamic turbulence, kinetic helicity is an ideal 
invariant and is known to cascade to smaller scales ll32ll33ll . 
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FIG. 1 1 : (color online) (a) Magnetic energy (upper curves) and ki- 
netic energy (lower blue curves) as a function of time, and (b) mag- 
netic helicity as a function of time. Solid lines correspond to DNS, 
dashed lines to 128'' a-model simulations, and dotted lines to 64"^ 
a-model simulations. 



Fig. [2]shows the spectrum of Hk during the hydrodynamic 
simulation. As with the energy, the alpha-model is able to 
capture the evolution of kinetic helicity in Fourier space up 
to fc ~ a~^. It seems that a Kolmogorov spectrum results 
for helicity iQEH, which implies that the relative helicity is 
weaker at small scales. 

The early stages of the growth of the magnetic field are in 
the "kinematic dynamo" parameter regime, involving expo- 
nential growth of the magnetic energy. This is shown in Fig. 
[Ts] which exhibits both the kinetic and magnetic energy as 
functions of time for the three runs. Though the three ener- 
gies as functions of time are offset, it is clear that the linear 
growth rates are close. At about t = 30, there is a saturation, 
close to a state in which on the average the energy is equipar- 
titioned approximately between kinetic and magnetic. After 
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FIG. 12: Spectrum of magnetic helicity for different times, for t — 
30, 35, 40, 45, 55, and 72.5; (a) DNS, (b) 128^ a-model, and (c) 
64'' a-model. The vertical lines indicate q^^. Note the cascade of 
magnetic helicity to large scales as time evolves. 



that, there are no significant variations in the evolution of the 
total kinetic and magnetic energy. 

Figs. ll6land[T71a.b show the negative of the magnetic he- 
licity, the mean square vector potential, and the mean square 
current density as functions of time. Though the agreements 
are not sharp, it is clear that the saturation levels and the times 
of saturation are both well approximated. Note that, in ac- 
cord with expectations l36il . the magnetic helicity acquires a 
negative value, opposite to the sign of the injected mechanical 
helicity. Note also that growth rates of both small scales (rep- 
resented by the square current) and large scales (represented 
by the square vector potential) are well approximated by the 
alpha model during the kinematic regime. 

While at t « 30 saturation in the exponential growth of 
magnetic energy takes place, the large scale modes continue 
growing, and at the end the magnetic field is dominated by 
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FIG. 13: Probability distribution function of relative kinetic belie- ^1^. 15: (color online) Kinetic energy (upper blue curves) and mag- 
ity. Solid lines correspond to DNS, dashed lines to 128^ Q-model netic energy (lower curves) as a function of time. Labels are as m 



simulations, and dotted lines to 64 Q-model simulations. 



Fig. EH 




FIG. 14: Spectrum of kinetic helicity. The Kolmogorov's slope is 
shown as a reference. The vertical lines indicate q^^. Labels are as 
in Fig. [HI 




FIG. 16: Negative of the magnetic helicity as a function of time. 
Labels are as in Fig. 1131 



large scales. While the mean square current density is con- 
stant after t =- 30, the squared vector potential keeps growing 
slowly. This behavior is even clearer in the evolution of the 
magnetic energy spectrum. 

Figs. ^]a,b show the evolution of the kinetic and magnetic 
spectra. The thick lines indicate kinetic spectra and the thin 
Unes the magnetic spectra. The vertical lines indicate the lo- 
cations of the two values of oT^ . In Fig. ^]a, the upper traces 
are the DNS spectrum at the time the seed field begins to grow, 
both for DNS and a-model simulations. The lower traces in 
Fig. ^]a show the magnetic spectrum at an early stage of 
its evolution. During this stage, the magnetic energy spec- 
trum peaks at small scales, and the a-model correctly captures 
the overall shape of the spectrum as well as the scale where 



the magnetic energy peaks. In the kinematic regime, all the 
magnetic /c-shells in Fourier space (up to /c < 12) grow with 
the same rate, and this feature is also well captured by the 
a-model simulations (not shown). This evolution is charac- 
teristic of small scale dynamos, as well as a k^l"^ slope in the 
magnetic energy spectrum at early times lls^ IstIi . Fig. 1181 b 
shows the late-time spectra, when approximate equipartition 
has been achieved. Note that as a result of helical dynamo ac- 
tion, a magnetic field at large scales (fc = 1) is generated (see 
Fig. [^b). The amplitude of this mode is in good agreement 
for both DNS and a-model simulations. 

Figure ^] shows surfaces of constant magnetic energy at 
t — 60, when the nonlinear saturation has already taken place 
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FIG. 17: (a) Mean square vector potential, and (b) mean square cur- 
rent density as functions of time. Labels are as in Fig. ^| 



but the large scale magnetic field is still growing. Thin and 
elongated structures can be identified in the magnetic field 
growing in the DNS. However, note that while these structures 
are present both in the DNS and in the alpha-model, in the lat- 
ter case the structures are thicker. This change is related to 
the filtering length a in the alpha-model. Similar results have 
been found in vorticity structures observed in simulations of 
hydrodynamic turbulence using the alpha-model fl. 

We thus conclude that there are few surprises in the dynamo 
simulations, at least for these values of rijv, and no glaring 
departures of the a-model predictions from the DNS results. 
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FIG. 18: (color online) Kinetic (thick blue lines), and magnetic en- 
ergy spectra (thin lines), at (a) f = 6 and (b) t = 70. Kolmogorov's 
hr^l'^ and Kazantsev's k?'^'^ spectra are shown as a reference. The 
vertical dotted and dashed lines correspond to the scales oT^ — 10 
and 20 respectively. Labels are as in Fig. 1131 



VII. SUMMARY; DISCUSSION 

Within the framework of rectangular periodic boundary 
conditions, we have examined four familiar three-dimensional 
MHD turbulence effects via the a-model and DNS. In every 
case, the principal large-scale features of these phenomena 
have been achieved with the a-model to acceptable accuracy. 
The savings in computer time achieved by the a-model runs 
have ranged from 256^/1283 = 8 to 256^/643 = 64, with- 
out considering extra saving in the time step from the CFL 
condition as the resolution is decreased. In no case has the 
a-model yielded results at significant variance with the DNS 
runs, which have been regarded as accurate. 

Other features of the DNS runs, such as the probability dis- 
tribution functions of the fluctuating quantities (such as local 
energy dissipation rates), have also been reproduced by the a- 



model as they were in two dimensions fl4ll . but we have not 
shown those results here because they are so similar to what 
was found in two dimensions. 

In Ref. fl4ll also the errors of the a-model computations 
were compared against under-resolved DNS. The behavior of 
the a-model in three dimensions is comparable to our previ- 
ous results, and therefore we refer the reader to our previous 
work for a detailed discussion about this topic. 

In conclusion, the MHD a-model can be considered to be 
validated, at least for the behavior of long-wavelength spectra 
in periodic boundary conditions. Its implementation in the 
presence of material boundaries stands as a next forbidding 
challenge. 
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FIG. 19: (color online) Surfaces of constant magnetic energy att = 
60 at 50% of its maximum value, for the DNS (above), and the 64"^ 
alpha-model (below). 
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